Probing the pairing symmetry and pair charge stiffness of doped t — J ladders 
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We perform the numerical equivalent of a phase sensitive experiment on doped t — J ladders. We apply 
proximity effect fields with different complex phases at both ends of an open system and we study the transport 
of Cooper pairs. Measuring the response of the system and the induced Josephson current, Density Matrix 
Renormalization Group calculations show how, depending on the doping fraction, the rung-leg parity of the pair 
field changes from minus to plus as the density of holes is increased. We also study the pair charge stiffness, 
and we observe a supression of the superconductivity in the region where static stripes appear. We compare our 
results with predictions from bosonization and renormalization group analysis. 

PACS numbers: 74.45.+c,74.50.+r,71.10.Pm 
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The 2-leg t — J ladder provides an example of how an ap- 
parently simple model of a strongly correlated many electron 
system can exhibit a rich variety of phenomena. UJJ This is il- 
lustrated by the schematic phase diagram reproduced in 
Fig.H] Here J/t is the ratio of the near neighbor exchange J 
to the near neighbor hopping t, and (n) is the site filling. For 
(n) = 1, one has a 2-leg Heisenberg ladder which has a spin 
gap A s w J/2. When holes are added, the system enters a 
spin gapped 'd-wave-4/ci?CDW' Luther-Emery phase for 
physical values of J/t. In this region there are power law pair- 
field and 4fcp-CDW correlations. For (n) > 0.5, the pairfield 
correlations are said to be 'd-wave' like because the correla- 
tion function of the singlet rung and singlet near neighbor leg 
pairfields is negative. For larger, unphysical values of J/t, 
the system phase separates. Density Matrix Renormalization 
Group calculations^ also provide evidence of a commensu- 
rate CDW phase at smaller values of J/t for (n) = 0.75 and 
0.5. At larger doping, where (n) < 0.5, there is a gapless Lut- 
tinger liquid phase and at low doping and large J/t values, an 
electron pairing region characterized by 's-wave'-like positive 
rung-leg pairfield correlations. 

Here we investigate the pairfield response of the 2-leg t — 
J ladder using a recently developed numerical technique[5] 
in which pairfields are applied to the ends of the ladder, as 
illustrated in Fig. |2j>). The Hamiltonian for the ladder has the 
usual form 
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Here, the operator c\ X(J creates an electron on rung i and leg 
A = 1, 2 with spin a, m\cr is the electron number operator 
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FIG. 1 : Schematic phase diagram of the t — J ladder as a function of 
J/t and site density (n), from Refs.^H)]. 



and Sia = c\ x ^Ci\. The Hilbert space excludes all states 
with double occupied sites. 

The coupling to the external pairfields shown in Fig. [2] is 
given by 

Hi = Ai (pt + /i.e.) + A 2 U+Pl + h.c)) , (2) 



where 
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creates a singlet pair on the first rung on the left end of the 
ladder and 

= ( c l-i,it c l,U — c i-i,i|. c i,iT) ^ 

creates a singlet pair on the last leg section the right end 
of the ladder (see Fig. |Zb). In the following we will set 
Ai=A2=<=l and use <j> to control the Josephson pair 
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FIG. 2: a) Comer configuration used in SQUID phase sensitive ex- 
periments, from Ref.(g|; b) Ladder geometry used in our calcula- 
tions; c) and d) Typical current response expected from d-wave and 
s-wave like superconducting orders, respectively. The sawtooth pro- 
file corresponds to perfect Andreev reflection, while the sine wave re- 
sponse describes the behavior as Andreev reflection becomes small. 



current through the ladder. The geometry of the external pair- 
field connections is similar to that used by Van Harlingen in 
his study of cuprate superconductors, jg] reproduced here in 
Fig.|2h), and to the arrangements proposed by Sigrist and Rice 
Q2D In the 'd-wave' like phase, we expect that there will be a 
7r phase shift associated with the sign difference between the 
rung and leg pairfields, as shown in Fig. [2] As we will discuss, 
the sawtooth curves in Figs[2j;) and d) represent the perfect 
Andreev reflection boundary condition limit and the Joseph- 
son sine waves represent the response in the case in which the 
boundary conditions evolve to zero Andreev reflection. ||5|] 

The Density Matrix Renormalization Group (DMRG) 
method ]8[] was used to calculate the ground state of the 
Hamiltonian H + H\. The term (O does not conserve the par- 
ticle number, but it preserves the total spin, and charge even- 
odd parity (number of fermions modulo 2). In order to work 
at a fixed density, the chemical potential fi has to be carefully 
calculated for each system size. Working in the grand canon- 
ical ensemble, as well as using complex numbers due to the 
Josephson phases, clearly limits our ability to go to very large 
systems. We have typically studied ladders of up to L = 24 
rungs, keeping a total 600 DMRG states, with a truncation 
error or the order of 10 -6 , or smaller. 

Away from the contacts, the current distribution is domi- 
nantly along the legs 
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and is independent of of the position of the rung i. Fig. |3^) 
shows DMRG results for L(J X ) for J/t = 1.0 and (n) = 
0.875 for ladders of various lengths L. As L increases, the 
current versus the phase <f> evolves into a 7r-phase shifted saw- 
tooth curve. This is what one would expect for a 'd-wave' like 
superconductivity state. That is, just as discussed in Ref.|@] 
for a one-dimensional electron liquid, as the length of the lad- 
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FIG. 3: Josephson current through t — J ladders of different lengths 
L, at a fixed density (n) = 0.875. The upper panel shows an in- 
creasing response which approaches a 7r-shifted d-wave sawtooth as 
L increases for J/t = 1.0. In the lower panel, for J/t = 0.5, the 
response approaches a Josephson sine wave form and decreases as L 
increases. 



der increases the boundary condition flows to one with perfect 
Andreev reflection when the pairing correlations are domi- 
nant. We previously found this in similar DMRG calculations 
for the one-dimensional Hubbard model [5] with an attractive 
interaction. In the present 2-leg t — J ladder, the difference is 
that the sawtooth is shifted by tt due to the d-wave nature of 
the pairing. 

In Ref.[9] bosonization and renormalization group methods 
were used to predict the scaling of the current with length, for 
a single leg Hubbard or t — J model. These methods can 
be readily extended to the two-leg case, based on the well- 
understood behavior of that model as studied, for example, in 
Ref . H . The rung boundary term, at the left-hand side of the 
system adds a term: 
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to the bosonized Hamiltonian. Here C is a constant (which 
can be taken to be positive) and <p(x) is the field arising from 
bosonization associated with the total charge degrees of free- 
dom which was labelled <j> +p in Ref. [2]. Taking into account 
the (normal) free end boundary condition (b.c), 9(0) = con- 
stant (where 8 is the field dual to ip), we see that this interac- 
tion has a renormalization group scaling dimension of 1/ (2K) 
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where K is the Luttinger parameter (referred to as K +p in 
Ref.|2[]). This is relevant for K > 1/2 and irrelevant for 
K < 1/2. In the relevant case, it is plausible that the cou- 
pling constant A renormalizes to infinity, changing the b.c. to 
ip(0) = 0, which corresponds to a perfect Andreev reflection 
b.c. We can test the consistency of this hypothesis by con- 
sidering the lowest dimension boundary operator which could 
appear in the effective fixed point Hamiltonian. For the single 
leg case this is a normal reflection term h.c. However, 

for the 2-leg ladder in the 'd-wave-4/ci?CDW phase studied 
in Ref. [2] this operator has exponentially decaying correlation 
functions and is strongly irrelevant. The most relevant bound- 
ary operator is (i/j^i/jr) 2 + h.c. oc cos(2y / 7r6') corresponding 
to normal pair reflection. This has a scaling dimension 2K 
and thus is relevant for K < 1/2 and irrelevant for K > 1/2. 
This strongly suggests that the renormalization group flow is 
from normal to Andreev fixed points for K > 1/2 and from 
Andreev to normal for K < 1/2, just like in the single leg 
case. 

The Josephson current can now be calculated by scaling ar- 
guments using these results. For K < 1/2 the pairing interac- 
tions scale towards zero as A(L) oc l/L 1 '^ 21 ^' 1 , so we ex- 
pect the current to scale to zero as (1/L)A(L) 2 oc l/L 1 ^^ 1 
and to have an approximately sinusoidal form at large L. For 
K > 1/2, as L — > oo, we expect to get a sawtooth current: 



J x L = 2vK{-l + (j)/-K), (0 < <j>< 2tt) 



(7) 



where v is the velocity of the gapless low energy charge ex- 
citations (referred to as v+ p in Ref.01). Here we have in- 
cluded the 7r-phase shift which arises from the fact that the 
leg pairing operator at the right end of the system leads to a 
term like H\ in the bosonized Hamiltonian but of opposite 
sign (and reduced magnitude). Eq. ^ can be derived using 
J = —2dE/d<f> and calculating E(<fi) using the low energy 
effective Hamiltonian of Eq. (2.6) of Ref.yl] with boundary 
conditions (^(0) = 0,cp(L) = \/rr— 4>/y/n. The corrections to 
this sawtooth behavior scale to zero with L and are controlled, 
at large L, by the leading irrelevant operator of dimension 2K. 
For instance, the critical current should behave as: 



J C L = 2vK[l - V(L) f ] 



where 



V(L) oc 1/L 



2K-1 
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is the renormalized normal pair reflection amplitude at scale 
L. The exponent / was argued to be / = 2/3 in Ref.l9(] by 
consideration of the non-interacting case, but it might be ex- 
pected to have a different value for the 2-leg ladder where the 
leading irrelevant operator corresponds to normal pair reflec- 
tion rather than normal electron reflection. We calculated K 
and v by studying the finite size scaling of the compressibil- 
ity and the excitation gap for S z = 0, for system sizes up to 
L = 56, and using 1500 DMRG states.[2] For J/t = 1 and 
(n) = 0.875, we find v = 0.313; K = 0.766; 2vK = 0.480, 
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FIG. 4: losephson current through t — J ladders of different lengths 
L, at a fixed density (n) = 0.25. The upper panel shows a sawtooth 
response for J/t = 3.0. In the lower panel, for J/t = 0.5, the 
response decreases with increasing system size, resembling a typical 
sine-like losephson shape. 



which is in good agreement with the limiting behavior of 
L(J X ) seen in Fig.fS^), that shows a slope of 0.484. 

From the phase diagram of Fig.Q] one sees that the parame- 
ter K decreases as J/ 1 is reduced. So, keeping the filling fixed 
at (n) = 0.875, we studied the behavior of L(J X ) for smaller 
values of J/t. Fig. [3}?) shows L(J X ) versus <\> for J/t = 0.5 
for ladders of various lengths. In this case, we find that the re- 
sponse maintains its n -phase shifted d-wave like character, but 
as L increases it evolves to a Josephson like sin (tp + ir) form. 
In addition, the amplitude decreases. This is what one would 
expect if the boundary conditions associated with the pairfield 
coupling evolve to zero Andreev reflection as L increases. 

Based upon the phase diagram of Fig.Q] we would have ex- 
pected that one would have had to go to a lower value of J/ 1 
for this to happen. In Ref.[2], the value of K was estimated 
by two different methods, obtaining K = 0.604 and 0.633 for 
J/t = 0.5 and (n) = 0.875. This is in apparent contradiction 
with the renormalization group arguments given above and the 
behavior seen in Fig.fSJ)). Thus, it may be that, for these pa- 
rameters where K is rather close to its critical value of 1/2, 
larger system sizes must be studied to see the asymptotic be- 
havior. Another possible explanation would be a non-trivial 
phase that exists for intermediate values of A, for a range of 
K. 
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FIG. 5: Josephson current through a t — J ladder with L — 16, 
for different values of J/t. D is defined as the slope of the current 
D — L\ ( J x ) |/0.l7r, for (f> = 0.97T. We show results as a function of 
the hole concentration 8 = 1 — (n) in the d-wave region S < 0.5. 
All quantities are in units where the hopping t = 1. 

In contrast to the 7r-phase shifted 'd-wave' response, a 
similar set of calculations for the case (n) = 0.25, shown 
in Figs|4^ and b, clearly exhibit an 's-wave' response. For 
J/t = 3.0, the L(J X ) curves exhibit the expected sawtooth 
form although the response appears essentially independent 
of system length L. For J/t — 0.5, when the CDW corre- 
lations are dominant, the flow is as expected to a Josephson 
sin <j> form associated with total normal reflection boundary 
conditions. While this is an unphysical region of the t — J 
ladder, it provides a further illustration of the way in which 
pairfield couplings can be used to explore the superconduct- 
ing response of lattice models. 

Finally, for a finite ladder, it is interesting to study the be- 
havior of L( J x )/4> versus the doping 6=1— (n) for various 
values of J/t. As seen from the schematic phase diagram 
shown in Fig. [T] at small values of J/t, one expects to cut 
through a CDW transition when S = 0.25. We define the ef- 
fective pair charge stiffness D as the slope of the Josephson 
current near = tt. Results showing D = L(J x )/(<fi — tt) 
versus 6 for several values of J/t are plotted in Fig. [5] As 
expected for small values of J/t, there is a sharp dip in D 
as 5 passes through 0.25. This dip decreases as J/t in- 
creases and one passes further away from the endpoint of 



the CDW line-phase. There are also dips associated with 
shifts in the doping AS = 1/16 corresponding to the addition 
of pairs. Smooth boundary conditions tend to remove these 
sharp features. II 1 Ofl Although one would need a detailed scal- 
ing study of LJ x /(<j) — tt) to determine the superfluid weight, 
the dip at 6 = 0.25 for J/t = 0.3 clearly shows that the t — J 
ladder exhibits a suppression in its superconductivity in the re- 
gion near where static stripes appear. Another aspect of this, 
which we have observed, is that when the boundary pairfield 
is relevant, the Friedel density oscillations generated by the 
open ends are reduced in amplitude. 

In summary, by applying pairfields we have carried out 
the numerical equivalent of a phase sensitive experiment 
on t — J ladders and studied the pair transport properties 
of a superconductor-Luther-Emery-superconductor (S-LE-S) 
junction. We have seen how the shape of the pair-current- 
phase characteristic can be used to determine the rung-leg par- 
ity of the superconducting pairs and its strength. One can also 
probe the effect of other phases on the superconducting re- 
sponse. 
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